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Thirteenth  Moving  Base 
Gravity /Gradlometer  Conference 
United  States  Air  Force  Academy 
Colorado  Springs,  Colorado 


Agenda 


Tuesday  12  February 

0730  -  Buses  depart  USAFA  Officers’  Club  for  Fairchild  Hall 

0745  -  Late  Registration  -  3rd  floor  Fairchild  Hall,  South  End 

0815  -  Wei come /Introduction  -  ILt  Donna  Warner 

0820  -  Opening  Remarks  -  Dr.  Charles  Martin  (DMA) 

0830  -  Gravity  Gradlometer  Survey  System  Program  Status  Review  - 
Mr.  Richard  Borgeson  (AFGL) 

0850  -  Review  of  Moving  Base  Gravity  Survey  System  (Navy)  at  Bell 
Aerospace  -  Mr.  Ernest  Metzger  (Bell) 

0920  -  Review  of  Moving  Base  Gravity  Gradlometer  Survey  System  (Air 
Force)  -  Mr.  Lou  Pfohl  (Bell) 

0950  -  Break 

1010  -  GGSS  Data  Analysis  Review  -  Dr.  Christopher  Jekell  (AFGL) 

1025  -  GGSS  Test  Planning  and  Test  Data  Treatment  -  Dr.  Warren  Heller 
(TASC) 

1055  -  Contribution  of  Terrain  In  Phase  II  Test  Area  to  Disturbance 

of  Gradients  Measured  Aloft  -  Mr.  C.  Lawrence  Bradley  (Geospace) 

1110  -  The  Template  Method  for  GGSS  Data  Reduction  -  Dr.  Jake  Goldstein 
(TASC) 

1145  -  Buses  Depart  Fairchild  Hall  for  USAFA  Officers'  Club 

1200  -  Luncheon  -  Officers'  Club 

1300  -  Buses  Depart  Officers'  Club  for  Fairchild  Hall 

1315  -  Synthetic  Data  Generation  for  Analysis  of  GGSS  Stage  II  Reduction 
Techniques  -  Dr.  Peter  Uglnclus  (NSWC) 

1345  -  Post  Mission  Adjustment  of  Gravity  Gradlometer  Data  - 
Mr.  G.  Scott  Peacock  (Geospace) 

1415  -  Frequency  Domain  Processing  of  Airborne  Gravity  Gradlometer 

Survey  Data  for  Derivation  of  the  Gravity  Disturbance  Spectrum  - 
Dr.  John  Hutcheson  (Bell) 
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1445  -  Break 
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1505  -  (a)  On  the  Problem  of  Astrogeodetic-Gradlometrlc  Vertical 
Deflection  Determination  in  Mountainous  Terrain 
(b)  A  High  Accuracy  Astrogeodetlc-Inertial  Data  Base  for  Tests 
of  the  Bell  Gravity  Gradlometer  -  Dr.  H.  Baussus  von  Luetzow 
(Army  -  ETL) 

1535  -  Analysis  of  the  High  Frequency  Spectrum  of  the  Anomalous 

Potential  -  Mr.  Anthony  A.  Vasslllou  (University  of  Calgary) 

1615  -  Buses  depart  Fairchild  Hall  for  Officers'  Club 

1630  -  Reception  -  USAFA  Officers'  Club 


Wednesday  13  February 

0745  -  Buses  Depart  USAFA  Officers'  Club  for  Fairchild  Hall 

0800  -  Two  Dimensional  Karhunen-Loeve  Smoothing  of  Gravity 
Gradlometer  Surveys  -  Dr.  Sam  Bose  (Applied  Science 
Analytics  Inc.) 

0830  -  Deterministic  and  Stocastic  Representation  of  Gravity  Gradient 
Using  Spherical  Harmonics  -  Mr.  Marvin  May  (NADC) 

0850  -  Review  of  Superconducting  Accelerometer  and  Gravity 
Gradlometer  Research  at  the  University  of  Maryland  - 
Dr.  Ho  Jung  Palk  (University  of  Maryland) 

0930  -  Break 

Presentation  of  Classified  Papers.  Only  persons  with 
a  security  clearance  of  at  least  SECRET  will  be  allowed 
to  attend  these  presentations. 

0950  -  Basic  Concepts  in  Real-time  Processing  of  Gravity  Gradlometer 
Data  -  Mr.  William  Chairetakls  (Sperry) 

1020  -  Inferring  Gradlometer  Environmental  Sensitivity  From  Long  Term 
Test  Data  -  Mr.  Donald  Benson  (Dynamics  Research  Corporation) 

1040  -  Short  Break 

1050  -  Techniques  for  Using  Multiple  Distributions  in  GGSS  Stage  II 
Data  Reduction  -  Mr.  Alan  Rufty 

1135  -  DoD  Executive  Session 

1215  -  Buses  leave 
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MINUTES  OF  THE  THIRTEENTH  ANNUAL 
MOVING  BASE  GRAVITY  GRADIOMETER  CONFERENCE 


United  States  Air  Force  Academy 
Colorado  Springs,  Colorado 

12-13  February  1985 


Tuesday  Morning  Session 


Welcome /Introduction  -  ILt  Donna  Warner  (AFGL) 

As  coordinator  of  the  two-day  conference,  Lt  Warner  opened  the  meeting  on 
behalf  of  the  co-sponsorlng  agencies,  the  Air  Force  Geophysics  Laboratory 
(AFGL)  and  the  Defense  Mapping  Agency  (DMA).  She  then  Introduced  Dr.  C. 

Martin,  Director  of  the  Advanced  Technology  Division,  DMA. 

Opening  remarks  -  Dr.  Charles  Martin 

Dr.  Martin  welcomed  participants  and  attendees.  He  made  comments  pertaining 
to  the  past,  present,  and  future  of  gravity  gradlometry.  He  made  note  of  the 
Navy  Initiatives  which  generated  Interest  In  the  field  five  years  ago  and  how 
DMA  Is  currently  cooperating  with  them  and  the  Air  Force  In  pursuing  research 
In  surface  and  airborne  measurement  of  gravity  gradients.  Dr,  Martin  forsees 
a  need  for  accuracy  of  several  orders  of  magitude  greater,  especially 
In  conjunction  with  the  space  program.  He  concluded  by  Introducing  Dr.  Randy 
Smith,  DMA's  Chief  Scientist  In  Geodesy  and  Geophysics. 

AFGL/DMA  Program  Review  -  Mr.  Richard  Borgeson  (AFGL) 

Mr.  Borgeson,  program  manager  of  the  Gravity  Gradlometer  Survey  System, 
presented  a  review  of  what  transpired  in  1984  on  the  GGSS  Program;  the 
test  site  was  selected,  the  aircraft  was  selected  and  tall  number  151384 
was  accepted  in  January  1985.  He  then  presented  future  milestones  on  the  program: 
lab  testing,  aircraft  testing,  and  system  acceptance  and  delivery,  among  others. 

Review  of  Moving  Base  Gravity  Survey  System  (Navy)  at  Bell  Aerospace  - 
Mr.  Ernest  Metzger  (Bell) 

Mr.  Metzger  highlighted  the  major  areas  of  activity  on  the  Navy  Gravity  Survey 
System  (GSS)  Program.  He  described  the  system  which  contains  1  gravity  sensor 
platform  on  which  the  gravity  Instruments  are  mounted,  1  GSS  electronics 
cabinet,  and  1  memory  processor.  He  went  Into  detail  describing  the  micropro¬ 
cessor  functions,  the  orientation  of  the  GGIs,  and  the  configuration  of  the 
electronics  cabinet.  Mr.  Metzger  reviewed  the  key  milestones  since  the  1984 
review  stating  that  major  progress  was  made  In  three  areas  1)  GGI  slip  ring 
assembly  2)  GGI  shock  test  and  3)  error  coeffeclent  calibration. 

Review  of  Moving  Base  Gravity  Gradlometer  Survey  System  (Air  Force)  - 
Mr,  A1  Jlrcltano  (for  Mr.  Lou  Pfohl  -  Bell) 

Mr.  Jlrcltano  first  emphasized  the  survey  and  test  programs  for  each  phase. 

Phase  I  (land)  and  Phase  II  (air)  of  the  GGSS  Program.  He  Identified  the 
test  methods  and  mentioned  that  data  would  be  run  through  the  Rolm  MVIOOOO 
computer  to  decide  whether  or  not  to  repeat  a  test.  He  showed  the  flow  of 
analysis  for  the  Phase  II  data  reduction  using  the  MV  100000.  Mr.  Jlrcltano 
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also  discussed  the  GGl  performance  monitor. 

Questions: 

Mr.  M.  Molny  (Sperry)  -  Which  vehicle  environment  (land  or  air)  is  more 
severe? 

Mr.  Jlrcltano  -  Both  are  tough.  Depends  on  many  features  including  the 
operation  of  the  vehicle. 

Mr.  Molny  -  Is  it  your  intent  to  establish  precise  control  of  the  track  or  to 
determine  precise  posltons? 

Mr.  Jlrcltano  -  We're  not  surveying  the  track,  but  with  GPS  we  will  get 
accurate  vehicle  positions. 

Mr.  A.  Carlson  (DMAHTC/GST)  -  Why  do  you  have  four  closely  spaced  astro  points 
in  Phase  It  land  testing? 

Mr.  Jlrcltano  To  determine  short  term  variation  of  Instrument  performance. 
GGSS  Data  Analysis  Review  -  Dr.  Christopher  Jekeli  (AFGL) 

Dr.  Jekeli  reviewed  Stage  II  of  Phase  II  GGSS  data  processing.  It  entails  the 
reduction  of  gravity  gradients  measured  at  aircraft  altitude  to  estimates  of 
the  surface  gravity  disturbance  vector.  He  mentioned  alternatives  to  optimal 
estimation.  Current  efforts  related  to  Stage  II  data  processing  are  turning 
to  comprehensive  simulations  to  study  the  various  processing  algorithms. 

Several  Issues,  such  as  data  management,  have  yet  to  be  addressed  in  detail. 

Questions : 

Mr.  P.  Uglncius  (NSWC)  -  Is  the  PSD  (on  your  graph)  for  the  vertical  gravity 
disturbance? 

Dr.  Jekeli  -  Yes 

GGSS  Test  Planning  and  Test  Data  Treatment  -  Dr.  Warren  Heller  (TASC) 

Dr.  Heller  discussed  the  technical  factors  driving  the  design  of  the  airborne 
test  program,  test-based  validation  of  GGSS  performance,  and  data  treatment 
scenarios.  These  Issues  are  related  to  our  knowledge  of  the  gravity  field 
in  the  north  Texas  test  area.  He  drew  inferences  as  to  the  accuracy  which 
the  gradlometer  Is  most  likely  to  achieve  In  that  location.  Dr.  Heller 
offered  a  summary  perpectlve  of  the  near  future  of  moving  base  gravity 
gradlometry. 

Ques  tlons : 

Mr.  Molny  -  Is  the  test  directed  toward  instrument  noise  or  gravity  estimation? 
Why  not  use  all  the  Information? 

Dr.  Heller  -  The  test  Is  directed  toward  gravity  estimation.  Part  of  the 
test  Involves  controlling  variables  and  understanding  them  quantitatively. 

A  discussion  between  Decker,  Ifetzger,  Molny,  and  Heller  followed  concerning 
use  of  all  the  data  points. 
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Mr.  T.  Sims  (NSWC)  -  Can  we  get  gravity  to  within  1  mGal  anywhere  in  the 
inner  zone? 

Dr.  Heller  -  We  always  have  edge  effects. 

Mr.  Sims  -  In  reality  is  there  an  inner  zone? 

Dr.  Heller  -  For  comparison  with  truth  data,  I  guess  that's  a  way  of  thinking 
of  it. 

Dr.  J.  Katz  (SSPO)  -  Is  there  any  reason  for  choosing  the  Texas  test  area? 

Lt  Fundak  (AFGL)  -  There  is  little  topographic  signature,  high  frequency 
gravity  signature,  and  good  truth  data. 

Dr.  Katz  -  Is  there  any  mandate  that  you  survey  over  land? 

Lt  Fundak  -  Yes,  the  truth  data  is  sufficiently  accurate  over  land. 

Mr.  Carlson  -  Where  does  the  active  model  come  from? 

Dr.  Heller  -  Deep  ocean  trench  data  and  high  alpine  area  data.  This  represents 
the  upper  bound  of  gravity  variation  in  any  area. 

Comment  by  Lou  Decker  (DMAAC)  -  Base  line  model  is  then  not  a  world-wide 
average  but  a  typical  local  model. 

The  Template  Method  for  GGSS  Data  Reduction  -  Dr.  Jacob  Goldstein  (TASC) 

Dr.  Goldstein  presented  an  overview  of  the  template  method  for  GGSS  data 
reduction.  The  method  consists  of  first  averaging  the  airborne  gradlometer 
measurements  on  suitably  chosen  zones,  and  then  processing  these  averages 
using  optimal  estimation  techlques.  He  also  discussed  why  the  appropriate 
quantity  to  be  estimated  during  GGSS  testing  Is  the  residual  gravity  disturbance 
vector,  which  Is  defined  as  the  departure  of  the  point  gravity  disturbance 
from  its  local  mean.  Dr.  Goldstein  presented  the  results  of  applying  the 
algorithm  for  estimating  the  vertical  component  of  the  residual  gravity 
disturbance  vector.  The  simulated  data  were  obtained  from  a  model  tuned  to 
the  test  area.  The  results  of  the  simulation  constitute  an  end-to-end  check 
on  algorithm  consistency. 

Questions : 

Dr.  Uginclus  -  Why  do  you  get  better  results  by  excluding  data? 

Dr.  Goldstein  -  Forcing  average  gains  on  some  data  may  prove  detrimental  to 
the  total  solution.  It  is  better  to  exclude  some  data  than  use  It. 

Dr.  Bose  (Applied  Sciences  Analytics  Inc)  -  How  Is  the  noise  characterized? 

Dr.  Goldstein  -  We  assume  white  noise.  Red  noise  can  be  taken  out  or  reduced 
to  the  same  level  as  white  noise. 

Dr.  Bose  -  When  W  (width  of  averaging  zone)  goes  to  zero,  do  we  get  optimal 
estimation? 


Dr.  Goldstein  -  Yes. 
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Mr.  Zorn  (Dynamics  Research)  -  How  are  errors  in  self  gradient  compensations 
accounted  for? 

Or.  Goldstein  -  Bell's  error  budget  Includes  these  errors. 

Mr.  Brozena  (NRL)  -  On  repeat  track  error  model  determination,  how  close  do 
altitude  and  position  have  to  be  navigated? 

Dr.  Heller  -  Bell  has  an  autopilot.  Initial  position  can  be  re-established 
by  resetting  the  barometric  altimeter,  and  GPS  data  processing  has  to  take  Into 
account  day  to  day  barometric  surface.  It  will  be  difficult  to  stay  within 
200m  In  all  coordlnants  for  an  800km  track. 

Discussion  ensued. 


Tuesday  Afternoon  Session 

Synthetic  Data  Generation  for  Analysis  of  GGSS  Stage  II  Reduction  Techniques  - 
Dr.  Peter  Uginclus  (NSWC) 

Dr.  Uginclus  began  with  the  explanation  that  NSWC  has  been  tasked  by  DMA  to 
provide  the  gravity  gradlometry  community  with  synthetic  data.  These  data 
will  provide  a  common  basis  for  assessing  and  evaluating  different  computational 
methods  that  are  proposed  for  reducing  the  gradient  measurements  In  the  GGSS 
program.  NSWC  has  constructed  for  this  purpose  two  different  data  bases: 

1.  Stochastic  model  with  known  statistics, 

2.  Deterministic  (non-stochastlc)  features  with  superimposed  stochastic 
short-wavelength  signals. 

The  models  consist  of  point  masses  and  dipoles.  Gravity  gradients  are  computed 
on  a  plane  600m  above  the  earth's  surface  on  a  grid  specified  by  the  survey 
requirements.  "Truth  data",  gravity  disturbances,  are  computed  on  a  5  X  3  km 
grid  on  the  surface.  Dr.  Uginclus  went  on  to  describe  in  more  detail  each  of 
the  two  models. 

Comments : 

Dr.  Bowln  (Woods  Hole)  -  One  can  get  a  better  model  by  using  smaller  masses 
so  as  not  to  destroy  long  wavelength  structure. 

Dr.  Uginclus  -  I  agree. 

Dr.  Heller  (referring  to  an  earlier  statement  by  Dr.  Uginclus)  -  The  TASC 
data  reduction  algorithm  and  the  Bell  data  reduction  algorithm  are  not 
competing,  Just  Independent. 

Dr.  Uginclus  -  Only  one  method  may  be  used  at  a  time,  therefore  In  the 
future  they  may  be  competing. 

Mr.  Decker  -  We  will  use  the  best  algorithm  available. 


Questions: 

Dr.  Jekell  -  Is  there  a  measure  of  non-statlonarlty  and  non-lsotropy  for 
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your  model? 

Dr.  Uglnclus  -  Not  now.  One  can  display  along-track  PSDs  or  take  the  average 
of  all  along-track  PSDs. 

Dr.  Martin  -  Are  all  existing  algorithms  based  on  prior  knowledge  of  the 
gravity  model? 

Dr.  Uglnclus  -  Yes.  The  two  algorithms  mentioned  this  morning  are. 

Dr.  Jordan  (Ceospace)  -  It  Is  very  common  to  use  prior  statistics.  If  you 
don't  have  much  data  then  prior  statistics  have  a  big  Influence.  If  you  are 
swamped  by  data  then  prior  statistics  have  no  big  effect.  Therefore  lack  of 
priors  are  no  serious  problem  to  us. 

Discussion  concerning  statistics  ensued  between  Jordan,  DeBra,  Martin,  and  Heller. 

Post  Mission  Adjustment  of  Gravity  Gradlometer  Data  -  Dr.  Stan  Jordan  (for 
Mr.  Scott  Peacock,  Geospace  ) 

Dr.  Jordan  described  a  technique  for  post-mission  adjustment  of  gravity 
gradlometer  data  that  uses  principal  component  decomposition  of  the  gravity 
model  covariance  matrix.  This  particular  form  of  the  Singular  Value 
Decomposition  algorithm  leads  to  a  very  Informative  transformation  of  the 
standard  Least-Squares  Collocation  technique.  This  transformation  allows 
calculation  of  the  estimates  without  matrix  Inversion,  and  identifies 
redundancies  and  near  redundancies  In  the  observations.  He  also  presented 
results  from  the  application  of  this  new  technique  to  gravity  gradient  analysis 
as  well  as  suggestions  for  further  development. 

Questions: 

Dr.  DeBra  -  It  is  nice  to  have  several  methods  (for  data  reduction)  but  have 
any  comparisons  been  made? 

Dr.  Jordan  -  No,  we  are  waiting  for  models  from  NSWC. 

Dr.  Bose  -  This  a  special  case  of  Kalman  filtering.  Why  not  use  the  general 
case  for  unbiased  states? 

Dr.  Jordan  -  We  only  use  biased  states. 

Mr.  Rufty  (NSWC)  -  Is  least-squares  incorporated? 

Dr.  Jordan  -  Yes,  it  Is  minimum  error  variance. 

Discussion  ensued  between  Martin,  DeBra,  and  Jordan  concerning  weights  and 
data  reduction. 

Frequency  Domain  Processing  of  Airborne  Gravity  Gradlometer  Survey  Data  for 
Derivation  of  the  Gravity  Disturbance  Spectrum  -  Dr.  Hutcheson  (Bell) 

Dr.  Hutcheson's  presentation  was  divided  Into  two  sections.  He  presented 
part  I  and  Mr.  Andrew  Grierson  presented  part  II.  Dr.  Hutcheson  explained 
that  the  measurement  of  gravity  gradients  using  the  Bell  Gradlometer  In  a 
moving  ground  vehicle  or  low  altitude  aircraft  necessitates  extraction  of  a 
wide  bandwidth  gravity  signature  In  the  presence  of  Instrument  and  environmental 
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noise.  He  described  an  optimal  technique  for  accomplishing  this.  The 
technique  incorporates  apriori  knowledge  of  these  processes  which  generate 
the  signature  and  noise  but  Is  adaptive  to  particular  mission  and  Instrument 
characteristics. 

Questions: 

Dr.  De  Bra  -  Why  limit  the  area?  (320  km  vs  500  km  square  blocks) 

Dr.  Hutcheson  -  We  are  waiting  for  a  routine  to  run  on  the  Data  General. 

Dr.  Jekell  -  Are  those  estimates  for  the  total  vertical  deflection? 

Dr.  Hutcheson  -  No,  we  make  a  correction  for  the  bias  and  the  trend. 

Mr.  Vasslllou  (Unlv.  of  Calgary)  -  What  about  data  extrapolation  beyond  the 
survey  area? 

Dr.  Hutcheson  -  It  is  difficult  to  do  In  two  dimensions. 

Dr.  DeBra  -  Will  subsequent  320  km  blocks  be  overlapped? 


Dr.  Hutcheson  -  It  Is  a  huge  Issue.  We  are  not  concerned  with  that. 

Part  II  -  Mr.  Grierson  (Bell) 

Questions: 

Dr.  Heller  -  What  kind  of  bandwidth  will  the  demodulator  have? 

Mr.  Grierson  -  Depends  on  actual  data.  We  should  try  to  extend  bandwidth  out 
to  the  signal  that  you  have.  If  the  signal  crosses  adjacent  noise  peaks, 
then  determination  will  be  poor. 

Dr.  DeBra  -  What  Is  the  sensitivity  of  the  weights?  How  much  change  can  you 
allow  before  you  need  a  new  set  of  weights? 

Mr.  Grierson  -  I  don't  know.  We  need  to  do  examples  for  various  bandwldths. 

a.  On  the  Problem  of  Astrogeodetlc-Gradlometrlc  Vertical  Deflection 
Determination  In  Mountainous  Terrain 

b.  A  High  Accuracy  Astrogeodetlc-lnertlal  Data  Base  for  Teats  of  the  Bell 
Gravity  Gradlometer  -  Dr.  H.  Baussus  von  Luetzow  (Army  ETL) 

Dr.  von  Luetzow  discussed  topics  a  and  b  pointing  out  problems  to  be  solved 
and  other  areas  of  concern  to  researchers  in  the  field. 

Analysis  of  the  High  Frequency  Spectrum  of  the  Anomalous  Potential  - 
Mr.  Anthony  Vasslllou  (Unlv.  of  Calgary) 

Mr.  Vasslllou  emphasized  that  the  knowledge  of  the  high  frequency  spectrum  of 
the  anomalous  gravity  potential  Is  Important  for  the  development  of  data 
processing  techniques  In  airborne  gradlometry.  This  spectral  Information  can 
be  obtained  by  analyzing  dense  gravity  data  by  FFT  methods.  An  attempt  has 
been  made  to  use  available  data  files  In  Canada  for  this  purpose.  He  explained 
that  the  data  considered  consisted  of  a  rather  uniform  set  of  (5*  X  5')  mean 
anomalies  in  most  of  the  continental  areas  and  four  sets  In  specific  areas 
with  an  average  data  spacing  of  about  1  km.  The  analysis  showed  that  the 


gravity  spectrum  has  to  be  resolved  to  about  degree  and  order  7000,  i.e.  to  a 
point  spacing  of  about  3  km,  to  resolve  the  gravity  anomalies  to  the  one  mgal 
level.  This  means  that  a  similar  data  spacing  would  be  required  for  airborne 
gradlometry  In  order  to  get  an  equivalent  resolution.  The  derived  spectrum 
was  then  used  to  determine  suitable  degree  variance  models  and  covariance 
functions  for  airborne  gradlometry.  Results  of  mean  5'  X  5'  gravity  anomalies 
covering  the  non-mount alnous  areas  of  Canada  Indicate  that  the  degree  variances 
of  the  anomalous  potential  decay  like  “3.6  gg  -3  gg  Kaula's  rule 

Implies.  This  would  mean  that  a  second-order  Gauss-Markov  model  agrees  better 
with  the  results  from  this  analysis  than  a  third-order  model. 

Ques  t Ions : 

Mr.  Zorn  -  How  do  results  relate  to  gravity  gradients. 

Mr.  Vasslllou  -  One  can  get  analytical  expressions  for  gradients  In  all  models. 

Comment,  Dr.  Baubsus  von  Leutzow  (Army  ETL)  -  Drawbacks  with  the  second  order 
Gauss-Markov  model  can  be  overcome  with  Taylor  series  and  finite  differences. 

A  discussion  among  Baussus,  Bose,  DeBra,  Heller,  and  Hutcheson  followed 
concerning  data  point  spacing  along  the  si  rvey  track  at  other  than  grid 
intersect  Ions . 


Wednesday  Morning  Session 

Two-Dimensional  Karhunen-Loeve  Smoothing  of  Gravity  Gradlometer  Surveys  - 
Dr.  Sam  Bose  (Applied  Sciences  Analytics  Inc) 

Dr.  Bose  opened  with  a  review  of  the  usual  assumptions  made  when  processing 
gravity  gradients.  Moving  base  gravity  gradlometry  (MBGG)  offers  a  much 
high*  f  resolution  compared  to  conventional  gravity  surveys.  The  increased 
speed  and  resolution  possible  with  MBGG  comes  with  a  price  tag;  that  of 
dealing  with  enormous  volumes  of  data  In  the  five  Independent  gradient  tensor 
components  of  the  gravity  gradlometer.  Gravity  gradients  are  not  fundamentally 
different  from  gravity  measurements  yet  their  processing  by  means  of 
conventional  techniques  Is  both  Inappropriate  and  inadequate.  The  simplistic 
approach  of  Integrating  a  single  gradient  component  along  the  track  of  the 
moving  platform  to  derive  a  gravity  profile  Ignores  the  Information  contained 
In  the  cross-track  and  gradients.  Dr.  Bose  presented  a  two-dimensional  signal 
processing  algorithm  tc  process  all  the  gradlometer  Jata  simultaneously.  The 
measurement  data  can  be  collected  at  any  altitude  on  or  above  the  surface  of 
the  earth  In  a  planar  grldded  form.  The  method  differs  from  conventional 
least  squares  techniques  In  that  It  accomodates  statistical  models  for  both 
the  gravity  signal  and  measurement  errors.  The  gravity  signal  model  Is 
derived  from  the  physical  theory  of  geodesy  and  Is  based  on  the  concept  of 
modeling  the  unknown  mass  distribution  below  the  survey  region  as  multiple  two- 
dimensional  white  noise  layers  representing  the  vertical  derivative  of  the 
disturbance  potential  to  any  order.  Such  a  model  requires  no  restrictions  of 
causality,  statlonarlty,  or  Isotropy.  The  problem  of  simultaneous  smoothing 
of  all  the  gradient  measurements  from  all  survey  traverses  in  the  region  is 
solved  by  representing  the  disturbance  potential  solution  as  a  Karhunen-Loeve 
expansion.  Estimating  the  gravity  field  simply  reduces  to  estimating  the 
Karhunen-Loeve  coefficients  which  are  uncorrelated  and  converging.  The 
Implementation  of  the  algorithm  requires  no  matrix  Inversions  and  can  handle 


large  amounts  of  data  in  a  computationally  efficient  manner 


Questions: 

Mr.  Benson  (Dynamics  Research)  -  What  is  lost  by  truncating  the  series? 

Dr*  Bose  -  Truncation  is  done  by  setting  M  »  K  (Nyqulst  limit). 

Mr.  Benson  -  Does  the  method  have  to  be  modified  for  real  time  estimation? 

Dr.  Bose  -  This  is  not  a  real  time  estimation  technique. 

Deterministic  and  Stocastic  Representations  of  Gravity  Gradients  Using 
Spherical  Harmonics  -  Mr.  M.  May  (NADC) 

Mr.  May  presented  expressions  for  the  anomalous  gravity  gradient  maps  in 
terms  of  spherical  harmonics.  Comparisons  among  different  coefficient  sets 
are  quantified.  A  stochastic  gravity  field  spherical  model  is  fit  and 
extrapolated  to  the  data.  From  this  the  corresponding  parameters  of  the 
planar  field  model  are  derived. 

Questions : 

Dr.  Jekeli  -  Why  are  there  no  gradients  (on  your  map)  in  the  Atlantic? 

Mr.  May  -  The  contour  Interval  is  5  E. 

Dr.  Jekeli  -  Is  the  expansion  to  degree  180? 

Mr.  May  -  Yea 

Lt  Fundak  -  What  is  the  tactical  filter? 

Mr.  May  -  The  filter  used  in  submarines  for  real-time  applications. 

Mr.  Molny  -  What  is  the  sensitivity  of  the  tactical  filter? 

Mr.  May  -  It  is  sensitive  to  low  frequencies. 

Mr.  Zorn  -  Is  the  equivalence  between  MAG  and  STAG  exact  for  the  planar  case? 

Mr.  May  -  Yes. 

Review  of  Superconducting  Accelerometer  and  Gravity  Gradiometer  at  the 
University  o*^  Maryland  -  Dr.  Ho  Jung  Palk  (Univ.  of  Maryland) 

A  sensitivity  goal  of  3  X  10^  E  Hz“^/2  been  defined  for  the  Gravity 
Gradiometer  Mission  (GGM)  at  NASA.  Dr.  Palk  and  his  associates  are  constructing 
a  three  axis  superconducting  gradiometer,  which  will  be  interfaced  with  a  new 
six  axis  superconducting  accelerometer  for  a  flight  on  the  Space  Shuttle. 

The  immediate  sensitivity  goal  is  2  X  10  E  which  will  be  Increased  by 

another  order  of  magnitude  for  the  mission.  Experiences  obtained  with  the 
prototype  single-axis  instrument  as  well  as  Innovations  in  superconducting 
technology  have  been  Incorporated  in  the  design  of  the  advanced  model  of  the 
gradiometer.  The  six-axis  accelerometer  will  measure  three  components  of 
linear  acceleration  and  three  components  of  angular  acceleration  simultaneously 
by  monitoring  the  motions  of  a  single  magnetically  levitated  proof  mass.  The 
outputs  of  this  accelerometer  will  be  used  to  isolate  the  platform  from  vibration 
and  to  compensate  for  dynamically  induced  errors  from  the  gradient  measurement. 


Dr.  Falk  discussed  the  principles  and  designs  of  these  Instruments  and  their 
development  schedules.  He  also  discussed  possible  applications  of  such 
sensitive  Instruments  for  fundamental  physics  experiments,  terrestrial  gravity 
surveys,  and  Inertial  navigation  surveys. 

Questions: 

Dr.  Martin  -  How  realistic  Is  your  work  schedule,  especially  the  flight  In  1989? 

Dr.  Falk  -  If  there  are  no  restrictions  on  funding  and  hiring  -  yes  It  Is 
realistic.  But  currently,  funding  Is  restricted  and  personnel  Is  limited. 

Dr.  Martin  -  Are  there  technological  problems? 

Dr.  Falk  -  Depends  on  the  goals  for  the  gradlometer.  The  6  axis  accelerometer 
Is  new  and  higher  risk. 

Dr.  Uglnclus  -  What  Is  the  most  sensitive  test  for  the  Inverse  square  law? 

Dr.  Falk  -  Depends  on  the  range  of  resolution.  At  a  few  kilometers  we  don't 
have  good  data. 

Basic  Concepts  In  Real-Time  Processing  of  Gravity  Gradlometer  Data  - 
Dr.  William  Chalretakls  (Sperry) 

Dr.  Chalretakls  discussed  techniques  for  obtaining.  In  real-time,  anomalous 
gravity  gradients  In  the  local  level  (North,  East,  Down)  frame  from  a 
gradlometer  consisting  of  three  carouselled  Bell  Gravity  Gradlometer  Instruments. 
He  first  described  a  technique  for  estimation  and  direct  compensation  of  the 
effects  of  low  frequency  Instrument  errors  on  the  gradients.  He  then  discussed 
the  technique  of  pair-processing  which  allows  the  real-time  estimation  of  the 
complete  set  of  gradients  from  only  two  Instruments.  He  further  presented 
the  Integration  of  these  two  techniques  Into  a  real-time  algorithm;  this 
algorithm  first  obtains  three  sets  of  pair-gradients  free  of  low  frequency 
Instrument  errors  and  subsequently  weights  these  continuously  (according  to 
most  recent  estimates  of  the  quality  of  Instrument  outputs)  to  obtain  a  final 
estimate  of  the  local  level  gradients.  He  claims  this  process  provides  the 
best  real-time  measurements  possible  without  the  aid  of  Information  external 
to  the  system. 

Inferring  Gradlometer  Environmental  Sensitivity  From  Long  Term  Test  Data  - 
Mr.  Donald  Benson  (Dynamics  Research  Corporation) 

Since  Mr.  Benson  did  not  offer  a  written  abstract  and  his  paper  was  classified 
SECRET,  no  notes  were  taken  on  his  presentation. 

Techniques  for  Using  Multipole  Distributions  In  GGSS  Stage  II  Data  Processing  - 
Mr.  Alan  Rufty  (NSWC) 

Mr.  Rufty  presented  an  extension  of  his  previous  work  on  the  approximation  of 
the  Earth's  gravity  field  by  point  masses.  The  treatment  of  point  multipole 
distributions  Is  accomodated  within  the  point  mass  formalism  previously 
developed.  In  this  connection  he  also  presented  the  accuracy  of  global  combined 
point  mass/dlpole  nonlinear  least  squares  fits  developed  to  model  tesseral 
gravity  fields.  Mr.  Rufty  described  an  extension  of  the  formalism  that 
develops  a  candidate  technique  for  the  operational  data  reduction  approach 
needed  In  support  of  the  Gravity  Gradlometer  Survey  System. 
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Questions: 

Dr.  Bose  -  Will  there  be  conversion  problems? 

Mr.  Rufty  -  You  have  to  choose  spacing  carefully  to  be  consistent  with  data 
processing. 

Dr.  Goldstein  -  How  time  consuming  Is  the  Implementation? 

Mr.  Rufty  -  Fairly  operationally  efficient.  We  are  Inverting  small  matrices. 

Lt  Fundak  -  What  Is  the  difference  between  high  pass/low  pass  vs  bandpass 
filters? 

Mr.  Rufty  -  None. 

Closing  -  Dr.  Martin  called  upon  Dr.  Dan  DeBra  of  Stanford  University  to 
offer  closing  remarks.  Dr.  DeBra  stated  that  he  was  glad  to  see  such  low 
errors  In  the  gravity  gradlometer  Instrument  (GGI).  But  he  also  pointed  out 
that  we  shouldn't  become  overly  confident.  There  Is  a  need  for  cross  checking 
using  LaPlace  transforms  and  carouselllng.  He  suggested  that  If  there  Is 
another  generation  of  GGIs  that  another  sensor  be  added  to  give  a  higher 
confidence  level.  Dr.  Debra  also  noted  that  the  principle  discussion  at  this 
year's  conference  revolved  around  data  fitting.  He  considers  this  a  very 
relevant  Issue  and  regards  Dr.  Peter  Uglnclus'  work  as  most  Important  In 
the  field.  Dr.  DeBra  mentioned  that  he  is  concerned  with  what  he  considers 
the  over-emphasis  on  statistical  modeling.  He  pointed  out  that  we  shouldn't  be 
diverted  by  the  amount  of  Importance  placed  on  statistics.  He  emphasized  the 
need  for  multiple  methods  for  fitting  the  data  and  that  we  need  Independent 
checks  for  conslstancy.  Dr.  DeBra  remarked  that  models  are  good,  but  not 
guarranteed.  Dr.  DeBra  made  a  quick  reference  to  data  reduction  and  survey 
techniques.  He  posed  the  question  of  whether  or  not  we  need  to  repeat  the 
track  at  the  end  of  the  survey.  He  stated  that  this  might  be  a  good  way  to 
avoid  losing  a  chunk  of  data  on  a  day  to  day  basis.  Other  points  addressed 
were  1)  data  processing;  should  we  use  all  the  data  or  not?  2)  practical  aspects 
of  the  need  for  aircraft  requirements;  we  need  to  distinguish  between  the  control 
of  the  aircraft  and  control  of  the  sensors.  In  general  Dr.  DeBra  considered 
the  program  to  be  progressing  well  and  making  strides  In  the  right  directions. 
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USE  STAGE  I  PROCESSING  YIELDING  1  SEC  SAMPLES  OF  6  GGI  OUTPUTS  IN  GGI 


Bell  Aerospace 


GRAVITY  GRADIOMETER  SURVEY  SYSTEM 
PHASE  2  LAND  TEST  (1-30  APRIL  1986) 
TESTD 


PHASE  II  AIRBORNE 


Bell  Aerospace 


1/6 


Bell  Aerospace 


Bell  Aerospace 


6GSS  DATA  ANALYSIS  (REVIEW) 


4F67Z-  IZ  VEZ  Fr 


HELL  AEROSPACE 

TASC 

NSWC 

DATA  ANALYSIS  WORKING 
GROUP  (D-A.W.G. ) 


PRIMARY  RESPONSIBILITY 

INDEPENDENT  ANALYSIS 

GRAVITY  MODEL 


COORDINATION 


I/? 


STAGE  II  DATA  PROCESSING 


LAND  VEHICLE  TEST 


J  ^TJ.dy^T^d.)  ,  j  = 


x.y. 


AIRBORNE  TEST 


DOWNWARD  CONTINUATION 

LARGE  DATA  SET 
INCOMPLETE  DATA  SET 


BOUNDARY-VALUE  PROBLEM 
NUMERICAL  INSTABILITY 

DIFFICULT  DATA  PROCESSING 

FINITE  EXTENT  R  ALIASING  ERRORS 


OPTIMAL  ESTIMATION 
LEAST-SQUARES  COLLOCATION 

NON-,  SUB-OPTIMAL  ESTIMATION 
INTEGRAL  FORMULAS 
SPACE  DOMAIN  L-S-C-  APPROXIMATION 
FREQUENCY  DOMAIN  L-S-C-  APPROXIMATION 


RECIPROCAL 

DISTANCE 

(STAG) 


LOGIOCSPAT  FREO  tCY/MI) 

Gravity  Disturbance  (CD)  power  spectral  densities  for  GGSS  Phase  II  test  area 


THE  ANALYTIC  SCIENCES  CORPORATION 


SP-4423-16 


GRAVITY  GRADIOMETER 
(GGSS)  TEST  PLANNING 
AND  TEST  DATA  TREATMENT 


12-13  FEBRUARY  1985 


Prepared  for; 

THIRTEENTH  MOVING  BASE  GRAVITY  GRADIOMETER  REVIEW 
United  States  Air  Force  Academy 
Colorado 


Prepared  by: 
W.G.  HELLER 


THE  ANALYTIC  SCIENCES  CORPORATION 
One  Jacob  Way 

Reading,  Massachusetts  01867 


THE  ANALYTIC  SCIENCES  CORPORATION 


ABSTRACT 


GGSS  TEST  PLANNING  AND 
TEST  DATA  TREATMENT 


W.G.  HELLER 


TASC  IS  RESPONSIBLE  FOR  INDEPENDENT  EVALUATION 
OF  THE  GGSS,  TEST  AND  SURVEY  PLANNING.  DATA  REDUC¬ 
TION  ALGORITHMS,  AND  OTHER  ISSUES  RELATING  TO  GGSS 
DEVELOPMENT.  DISCUSSED  IN  THIS  PRESENTATION  ARE 
TECHNICAL  FACTORS  DRIVING  THE  DESIGN  OF  THE  AIRBORNE 
TEST  PROGRAM.  TEST-BASED  VALIDATION  OF  GGSS  PERFOR¬ 
MANCE  AND  DATA  COLLECTION  SCENARIOS.  THESE  ISSUES 
ARE  RELATED  TO  OUR  KNOWLEDGE  OF  THE  GRAVITY  FIELD  IN 
THE  NORTH  TEXAS  TEST  AREA  AND  INFERENCES  ARE  DRAWN 
ABOUT  LIKELY  ACCURACY  WHICH  THE  GRADIOMETER  WILL 
ACHIEVE  THERE.  A  SUMMARY  PERSPECTIVE  OF  THE  NEAR 
FUTURE  OF  MOVING  BASE  GRADIOMETRY  IS  OFFERED. 
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PROJECTED  GGSS  SURVEY  PERFORMANCE  USING 
DIFFERENT  GRAVITY  FIELD  MODELS 
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FOREWORD 


This  document  contains  materiai  used  in  a 
presentation  given  by  The  Analytic  Sciences 
Corporation.  The  material  is  not  intended  to  be 
self-explanatory,  but  rather  should  be  considered 
in  the  context  of  the  overall  presentation. 
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